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Abstract 



We consider one-dimensional chain of coupled linear and nonlinear oscillators 
> , 
^O ! with long-range power wise interaction defined by a term proportional to l/|n — 

m . 

'nJ" ■ TTT-I"^^. Continuous medium equation for this system can be obtained in the so- 



called infrared limit when the wave number tends to zero. We construct a transform 
operator that maps the system of large number of ordinary differential equations 
of motion of the particles into a partial differential equation with the Riesz frac- 
KJi \ tional derivative of order a, when < a < 2. Few models of coupled oscillators 

S . are considered and their synchronized states and localized structures are discussed 

in details. Particularly, we discuss some solutions of time-dependent fractional 
Ginzburg-Landau (or nonlinear Schrodinger) equation. 
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1 Introduction 

Although the fractional calculus is known for more than two hundred years and its developing 
is an active area of mathematics, appearance and use of it in physical literature is fairly recent 
and sometimes is considered as exotic. In fact, there are many different areas where fractional 
equations, i.e., equations with fractional integro-differentiation, describe real processes. Between 
the most related areas are chaotic dynamics [Ij , random walk in fractal space-time [2J and random 
processes of Levy type [3llll|5l|6]. The physical reasons for the appearance of fractional equations 
are intermittancy, dissipation, wave propagation in complex media, long memory, and others. 
An important area of the application of long-range interaction includes collective oscillation 
and synchronization in physics, chemistry, biology, and neuroscience actively studied recently 
[71 El El- Beginning from the pioneering contributions by Winfree |10] and Kuramoto jllj . 
studies of synchronization in populations of coupled oscillators becomes an active field of research 
in biology and chemistry. The synchronization can also be considered in complex oscillatory 
medium, where each site (element) performs self-sustained oscillations. A good physical and 
chemical example is the oscillatory Belousov-Zhabotinsky reaction [121 [T3l [TT] in a medium 
where different sites can oscillate with different periods and phases. 

Complex Ginzburg-Landau equation is canonical model for oscillatory systems with local 
coupling near Hopf bifurcation. Recently, Tanaka and Kuramoto [14] have shown how, in the 
vicinity of the bifurcation, the description of an array of nonlocally coupled oscillators can be 
reduced to the complex Ginzburg-Landau equation. In Ref. [15] . a model of population of 
diffusively coupled oscillators with limit-cycles is described by the complex Ginzburg-Landau 
equation with nonlocal interaction. 

Nonlocal coupling means interaction with oscillators or other objects distanced arbitrary 
far from each others. The corresponding interaction potential could be with a characteristic 
length Refs. [161 HTl [T5] or without it. The long-range interaction that decreases as l/lxl""*"^ 
with < a < 2 was considered in Refs. [THl dSl 1201 Ell 122] with respect to the system's 
thermodynamics and phase transition. It is also shown in p3] that using the Fourier transform 
and limit for the wave number A; — )• 0, the long-range term interaction leads under special 
conditions to the fractional dynamics. 

In the last decade it is found that many physical processes can be adequately described by 



equations that consist of derivatives of fractional order. In a fairly short period of time the list 
of such applications becomes long and the area of applications is broad. Even in a concise form, 
the applications include material sciences [231 [251 [111 [2Z| , chaotic dynamics [1] , quantum theory 
[251 [23 [30] 5 physical kinetics [H [31 [321 [31] , fluids and plasma physics [331 [33] j and many others 
physical topics related to wave propagation [35J , long-range dissipation [36] , anomalous diffusion 
and transport theory (see reviews [H [21 [231 [31 [3Z])- Particularly, fractional Ginzburg-Landau 
equation was suggested in [3S1 [33 130] • 

In this paper we would like to strength the state that the one-dimentional chain of parti- 
cles with long-range power type interaction can be asymptotically described by the continuous 
medium equation with fractional space derivative. We introduce, under some conditions, a corre- 
sponding transform operator that map the system of large number ordinary differential equations 
of particles motion into the equation with partial and fractional derivatives. It is worthwhile to 
mention that possibility of such an operation has an immediate physical consequence: fractional 
dynamics equation has solutions that describe coherent structures with power-like tails. 

In Section 2, we introduce the transform operator. In Section 3, some particular solutions are 
derived with a constant wave number for the fractional Ginzburg-Landau (FGL) equation that 
can be interpreted as synchronization in the oscillatory medium. In Section 4, we consider forced 
FGL equation for isochronous case. In Section 5, we derive solutions of the fractional Ginzburg- 
Landau equation near a limit cycle. These solutions are interpreted as coherent structures in 
the oscillatory medium with long-range interaction. 

2 Long-range interaction of oscillators 

Consider a system with Hamiltonian 



«= E 



11^°° 1 

l^ul+ -go Yl |„_^|l+« ""nUm + V{Un) 



2 ^-^ \n — m\ 

m=—oo 
ra^n 



(1) 



where g^ and a are some constants, and u„ = Un{t). The corresponding equations of motion are 

m=—oo ' ' 



with V'{u) = dV{u)/du. Let us introduce notations 

-, +00 
Hint = 1^90 ^ |n-m|"(^+°) n„Um, (3) 



2' 

m=—oo 

and 

+00 



f^un= Yl \n-m\-^^+^^Um. (4) 

m=— 00 

The standard transformation from the set of equations ([2]) to the continuous medium equations 
is in replacements 

Un{t) = u{Xn,t) ^ u{x,t), Xn = nAx, (5) 

where n is a positive integer number, Ax is a distance between particles, and the limit Ax — t- 
is applied. 

A specific form of Hint and the corresponding term T^ in the equation of motion ([2]) create 
a possibility to present the continuous medium equations in a form that consists of fractional 
derivatives. For this goal, we use the notation 

u{x, t) = — dk u{k, t) e^^^' = T-^{u{k, t)}, 

/ + 00 
dx Z{x, t) e"*''^' = J'{u(rc, t)}, (6) 

-oo 

where J-'{u{x,t)} is Fourier transform of u{x,t) with respect to x. Then, replacing ([5]) in ([2]), 
we obtain TaUn{t) — )• Tau{x,t), and after applying Fourier transform to ([2]), 

d'^u{k,t) 
where 



+ go[Ta{k)-Tamu{k,t)+ T{V'{u{x,t))} = 0, (7) 



T^ik) = 2T{-a) cos{7ra/2)\kr + 2 ^ ^^" + ^ ^""^ (-fc^)", |fc| < 1, a ^ 0, 1, 2,3..., (8) 

n=0 ^ '^'^' 

and Tq,(0) = 2^(1 + a) and C is the Riemann zeta- function. Function Ta{k) can also be presented 
in the form 

n=— CO ' ' n=l 



where Lia{z) is a polylogarithm function 

Li;3(e^)=r(l-/3)(-zf-i+f;^i^^z", |z|<27r, /3/ 1,2,3.... (10) 

n=0 

After substitution ([8]) into ([7]) obtain 



rn(fc, t)- 2,0 E^^^^ 



5i2 

where 

Oc, = 2r(-a)cos(7ra/2), a/ 0,1,2,.... (12) 

From Eq. (jlip we obtain the equation for the field u{x,t) using the definition ([6]) 

^^ = V'iuix, t)) + 5oa.^n(x, t) - 250 f; ^^^^^^^£!r-(-' *)' (« ^ 0, 1, 2...). 

(13) 
Here, we use the connection between Riesz fractional derivative and its Fourier transform 



o\x\°' o\x\'^ 

Such construction that leads to the equation with fractional derivative appears also in |23j . 
The properties of the Riesz derivative can be found in [431 [44l I45j . Its another expression is 

-f^u(x, t) = — i-— - (PXx, t) + P^n(x, t)) , (15) 

a|a;|° 2cos(7ra/2) 

where a 7^ 1,3,5..., and P^ are Riemann-Liouville left and right fractional derivatives 



r(n - q) dx"^ J_^ (x - ^)"-«+i ' 

" ^ ^ r(n - a) Sx'^ 7^ (C - 2;)"-"+i ^ ^ 

where n — 1 < a < n. 

The obtained expressions (|lip and (|13p can be simplified in the so-called infrared approxi- 
mation A; — )• 0. Then Eq. pip can be rewritten as 



—u{k,t)=T{V'{u{x,t))}+goTa{k)u{k,t), k ^ 0, (1< q < 3), (17) 

where we use 7^(/c) as an approximation expression for Ta{k): 

7;(/c) = 2^(-a)cos(7ra/2)|fc|"-C(a-l)A:^ (1< q < 2), (18) 



T2{k) = const InA;^, (a = 2), 
Ta{k) = Jok\ (a > 2). 
The expression for Ta{k) can be considered as a Fourier transform of an operator Ta{x): 

faix) = T-HTaik)} = aag^ + C(a - 1)^. (19) 



Applying inverse Fourier transform to ()17p . we also obtain 

■^u{x,t) = V'{u{x,t))+goTa{x)u{x,t), (l<a<2). (20) 

This equation can be considered as an equation for continuous oscillatory medium with 1 < a < 2 
in the infrared {k — )• 0) approximation. 
Let us note that (jlSp has a scale Iq: 

l^^ = \2T{-a) cos(TO/2)/C(a - l)|^/(2--a) (21) 

such that nontrivial expression Ta{k) ~ |^|" appears for k <^ k^. In the following, we consider 

the case k <^ kQ, i.e. 

7;(x)~a„-f— . (22) 

o\x\" 

The expression ([8]) for Ta{k) maps the system ([2]) for large number of oscillators into the con- 
tinuous medium equation ^ in the limit ([5]). The approximate expression Ta{k) (fTHj) for Ta{k) 
is used to construct an operator Ta{x) that stands instead of Ta in dU and that appears in the 
final approximate equation (j20p . We will call Taix) the transform operator. The reason for that 
is that it can be used for in more broad sense than just for the Hamiltonian ([1]). Here are few 
examples. 

The dynamical equations can be gain the term of the type ([3]) in complex dispersive media 
[26\ [35} I31j . The corresponding approximate equations with fractional derivatives appear for 
tracer dynamics in the presence of convective rolls |41j, and for the equation for surface wave 
interaction in ocean |42) . A specific type of equations is widely used for the oscillatory media: 

-7-Zn{t) = F{Zn) + go ^ Jain-m){Zn- Zm), (23) 

m=—oo 

where Zn is the position of the nth oscillator in the complex plane, and F is a force. As an 
example, for the oscillators with a limit cycle, F can be taken as 

F{z) = {l + ia)z-{l + ib)\z\'^z. (24) 
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The nonlocal interaction is given by the power function 

Ja{n) = |nr"~\ (25) 



Similarly to @. Using the notations (jlSp and (|19p . we can rewrite the approximate equation 
with fractional derivative in the form 

^^^ = F{Z{x,t))+gofe.{x)Z{x,t), (l<a<2), (26) 

where we replace Zn{t) — )• Z{x,t) and consider A; — )• 0. 

Using, for example F(z) = and expression ()22p for Taix), we arrive to the fractional kinetic 
equation: 

^Z(x, t) = goaa-gTT^Zix, t), (1< Q < 2) (27) 

that describes the fractional super diffusion [3l [3T| Hj. For F(z) in the form (|24p . we obtain 
fractional Ginzburg-Landau equation suggested in |381 [39l HO] that will be considered in the 
next section. 

Finally, consider 

I +°° 1 

Hint = ^90 Yl W;~^;^\T+^S{Un) giUm), (28) 

where g{u) is some function of u, instead of ([3]). The corresponding generalization of (|20p is 
easily obtained as 

^ ""dt/^ =V'{u{x,t)) + go g'{u{x,t))fa{x) g{u{x,t)), k ^ 0, (1< a < 2), (29) 

where g'{u) = dg{u)/du and Ta{x) is the same as in (fT9]) or ([22]) . 



3 Fractional Ginzburg-Landau (FGL) equation 

The one-dimensional lattice of weakly coupled nonlinear oscillators is described by 

— Z„(t) = (1 + ia)Zn - (1 + ib)\Zn\'^Zn + (Ci + iC2)(2:„+i - 2Zn + Zn-l), (30) 

at 
where we assume that all oscillators have the same parameters. A transition to the continuous 
medium assumes [8J that the difference Zn+i — Zn is of order Ax, and the interaction constants 
ci and C2 are large. Setting ci = 5f(Ax)~^, and C2 = (7c(Ax)~^, we get 

-Z = (1 + m)Z - (1 + i5)|Z|2z + 5(1 + ic)^Z, (31) 



which is a complex time-dependent Ginzburg-Landau equation. The simplest coherent structures 
for this equation are plane- wave solutions [8J: 

Z{x, t) = R{K)exp[iKx - iuj{K)t + 6*0], (32) 

where 

R{K) = {l-gK^f'\ u{K) = {b-a) + {c-h)gK\ (33) 

and ^0 is an arbitrary constant phase. These solutions exist for gK^ < 1. Solution (j32p can be 
interpreted as a synchronized state |8]. 

Let us come back to the equation for nonlinear oscillators (|23|) with F{z) in Eq. (|24|) and 
long-range coupling (j25|) : 



— Zn = (1 + ia)Zn - (1 + ib)\Zn\'^Zn + 50 V] 1 U;TTi^n " ^m), (34) 

at ^-^ \n — TTi "+^ 

ra^n 

where z„ = 2„(t) is the position of the nth oscillator in the complex plane, 1 < a < 2. The 
corresponding equation in the continuous limit and infrared approximation can be obtained in 
the same way as (pO]l : 

-Z = (1 + ia)Z - (1 + ih)\Z\^Z + g{\ + ^c)——Z, (35) 

where g'(l + ic) = gociai and 1 < a < 2. Eq. psp is a fractional generalization of complex time- 
dependent Ginzburg-Landau equation ([3T|) (compare to ([20]) 1. Here, this equation is derived in 
a specific approximation for the oscillatory medium. 
We seek a particular solution of (j35p in the form 

Z{x,t) = A{K,t)e'^''. (36) 



Eq. (|36p represents a particular solution of (|35p with a fixed wave number K. 
Substitution of ([36]) into ([35]) gives 

— ^(K, t) = (1 + ia)A - (1 + i6)M|2^ - g(l + ic)\K\'^A. (37) 

at 

Rewriting this equation in polar coordinates 

A{K,t) = R{K,t)e''^^^''\ (38) 



we obtain 



— = {l-g\KnR-R^, - = (a-cg\Kn-bR\ (39) 



The limit cycle here is a circle with the radius 

R = {1 - glKl")^/"^ , g\K\"<l. (40) 

Solution of (|39|) with arbitrary initial conditions R{K,0) = Rq, 9{K,0) = 9q is 

R{t) = Roil - g\Kn'^' [rI + il-g\Kr- i?2)e-2(i~5l^r)*) ~'^' , (41) 

9{t) = -^ In [(1 - g\Kn-' [R^ + (1 - g\K\- - i?g)e-2-*)] - u;^{K)t + 9o, (42) 



where 



ojaiK) = {b - a) + {c - b)g\K\" , 1 - g\K\'^ > 0. (43) 



This solution can be interpreted as a coherent structure in nonlinear oscillatory medium with 
long-range interaction. 

URl = l- glKl"", c/|i^|" < 1, then Eqs. gl]) and i^ give 

Rit) = Ro, 9{t) = -Ua(K)t + 9o. (44) 



Solution (|44p means that on the limit cycle (j40p the angle variable 9 rotates with a constant 
velocity Ua{K). As the result, we have the plane- wave solution 

Z{x, t) = (1 - ^|K|°)V2gii^x-i^<,(x)t+i0o ^ l-g\K\">0, (45) 

which can be interpreted as synchronized state of the oscillatory medium. 
The plane-wave solution (j45p is stable if parameters a, b, c and g satisfy 



< 1 - glKl" < a/b - {c/b)g\K\'' < 3(1 - ^1^1"). (46) 

Condition (I46p defines the region of parameters for plane waves where the synchronization exists. 
For initial amplitude that deviates from ()40p . i.e., i?Q 7^ 1 — ^li^l", an additional phase 
shift occurs due to the term which is proportional to b in (j42p . The oscillatory medium can be 
characterized by a single generalized phase variable. To define it, let us rewrite (|39p as 



±lnR = {l-g\Kn-R\ ±9 = {a - cg\Kn - bR\ (47) 



From (|47p . we obtain 

^<^ = -a;,(if). (48) 

where 

(p{R,e) = e-blnR (49) 

is the generahzed phase [8j, and uja{K) is defined by (03]). 

Eq. (|48|) means that generahzed phase 4>{R, 6) rotates uniformly with constant velocity. For 
g|-fC|" = (6 — a)/{h — c) < 1, we have the lines of constant generalized phase. On {R,9) plane 
these lines are logarithmic spirals 6 — b\nR = const. The decrease of a corresponds to the 
increase of K. For the case 6 = instead of spirals we have straight lines (j) = 6. 

Energy propagation can be characterized by the group velocity Va,g = duJa{K)/dK. From 
Eq. (03]), we obtain 



v^,g = a{c-h)g\K\^-\ (50) 



For \K\ < Ki = (a/2) ", we get \va,g\ > \v2,g\- The phase velocity is 

v^,ph = 0Ja{K)/K = {c-b)g\K\''-\ (51) 

For \K\ < K2 = 2"~^, we have |va,p/i| > |'V2,p?t|- Therefore long-range interaction decreasing as 
|2.|-(a+i) -^[iii 1 < a < 2 leads to increase the group and phase velocities for small wave numbers 
{K — )• 0). Note that the ratio Va,g/va,ph between the group and phase velocities of plane waves 
is equal to a. 

4 Forced FGL equation for isochronous case 

In this section, we consider FGL equation ()37|) forced by a constant E (the so-called forced 
isochronous case {b = 0) [8j): 

8 

—A = (l + ia)A-\A\'^A-g(l + ic)\KrA-iE, (lm^ = 0), (52) 

ot 

where A = A{K,t), and we put for simplicity 6 = 0, and K is a. fixed wave number. Our 

main goal will be transition to a synchronized states and its dependence on the order a of the 

long-range interaction. 

Numerical solution of Eq. (|52p was performed with parameters a = 1, g = 1, c = 70, 

E = 0.9, K = 0.1, for a within interval a G (1;2). The results are presented on Fig. 1. For 

10 



ao < a < 2, where oq ~ 1.51..., the only stable solution is a stable fixed point. This is region 
of perfect synchronization (phase locking), where the synchronous oscillations have a constant 
amplitude and a constant phase shift with respect to the external force. For a < oq the global 
attractor for (|52p is a limit cycle. Here, the motion of the forced system is quasiperiodic. For 
a = 2 there is a stable node. When a decreases, the stable mode transfers into a stable focus. 
At the transition point it loses stability, and a stable limit cycle appears. As the result, we have 
that the decrease of order a from 2 to 1 leads to the loss of synchronization (see Figure 1). For 
a = 2.00, and a = 1.60, we see that in the synchronization region all trajectories are attracted 
to a stable node. 

In Fig. 1, for a = 2.00, and a = 1.60, we see that in the synchronization region all trajectories 
are attracted to a stable node. For a = 1.54, a = 1.52, a = 1.50, we see that a stable limit 
cycle appears via the Hopf bifurcation. For a = 1.54, and a = 1.52, near the boundary of 
synchronization the fixed point is a focus. For a = 1.50, and a = 1.4, the amplitude of the limit 
cycle grows, and synchronization breaks down. 

The oscillator medium can be characterized by a single generalized phase variable (j49p . We 
can rewrite ()49p as 

(t>{X, Y) = arctan(y/X) - - ln(X2 + Y^), (53) 

where X and Y are real and imaginary parts of A{K, t). For E = 0, the phase rotates uniformly 

^J = -uJa{K)=a-gc\K\'', (54) 

at 

where uja{K) is gived by (j43p with 6 = 0, and can be considered as a frequency of natural 
oscillations. For E ^ 0, Eqs. (f^2|) and ([55]) give 

— (p = —uja{K) — E cos (j). (55) 



This equation has an integral of motion. The integral is 



I = 2|a;2 - ^2 1 -1/2 avctan(sgn{u; - E) 



io-E 



-1/2 



tan((/)(t)/2)) +t. (56) 



oj + E 

These expressions help to obtain the solution in form (j38p for forced case (j52p keeping the same 
notations as in (f55|) . For polar coordinates we get 

'^ = {l-9\KnR-R'-Esine, ^^ = ia-cg\Kn-^. (57) 
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Figure 1: Approaching to the bifurcation point a = ao = 1.51... and transformation to 
the hmit cycle of solution of forced FGL equation for isochronous case with fixed wave 
number K = 0.1 is represented by real X{K,t) and imaginary Y{K,t) parts of A{K,t). 
The plots for orders a = 2.00, a = 1.60, a ^1.54, a = 1.52, a = 1.50, a = 1.40. 



Numerical solution of (|57|) was performed with the same parameters as for Eq. (|52|) . i.e., a = 1, 
g = 1, c = 70, E = 0.9, K = 0.1, and a within interval a £ (1, 2). The results are presented in 
Figs. 2 and 3. 

The time evolition of phase 6(K,t) is given in Fig. 2 for a = 2.00, a = 1.50, a = 1.47, 
a = 1.44, a = 1.40, a = 1.30, a = 1.20, a = 1.10. The decrease of a from 2 to 1 leads to the 
oscillations of the phase 9{K,t) after the Hopf bifurcation at ao = 1.51..., then the amplitude 
of phase oscillation decreases and the velocity of phase rotations increases. 

The amplitude R{K,t) is shown on Fig. 3 for a = 1.6, a = 1.55, a = 1.55, a = 1.51, 
a = 1.50, a = 1.45, a = 1.2. The appearance of oscillations in the plots means the loss of 
synchronization. 

5 Space-structures from FGL equation 

In previous sections, we considered mainly time-evolution and "time-structures" as solutions for 
the FGL equation. Particularly, synchronization process was an example of the solution that 
converged to a time-coherent structure. Here, we focuse on the space structures for the solution 
of FGL equation (|35p with b = c = and the constants oi and a2 ahead of linear term 

^Z = {ai + ia2)Z-\Z\''Z + g^Z. (58) 

ot o\x\°' 

Let us seek a particular solution of (j58p in the form 

Z(x,t) = i?(x,t)e^^W, R*{x,t) = R{x,t), e*{t) = e{t). (59) 

Substitution of ([^^ into ([^5|) gives 

§.^B. = a,R-R^-S^R, !<,(*) =„,. (60) 

1/2 

Using 6{t) = a2t + ^(0), we arrive to the existence of a limit cycle with Rq = a^ . 

A particular solution of ()60p in the vicinity of the limit cycle can be found as an expansion 

R{x,t) = Ro + eRi + e^R2 + ..., (e <Cl). (61) 

Zero approximation Rq = a^ satisfies (|60p since 9°/9|x|"l = 0, and for Ri = Ri{x,t), we have 

|H, = -2a,iJ.+,^H,. (62) 
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Figure 2: Phase 9{K,t) for K = 0.1 and a = 2.00, a = 1.50, a = 1.47, a = 1.44, 
a = 1.40, a = 1.30, a = 1.20, a = 1.10. The decrease of order a corresponds to the 
clockwise rotation of curves. For upper curve a = 2. For the most vertical curve a = 1.1. 
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Figure 3: Amplitude R{K,t). The upper curve corresponds to a = 2 for all plots. The 
lower curves correspond to a = 1.6, a = 1.55, a = 1.51, a = 1.50, a = 1.45, a = 1.2. The 
appearance of oscillations on the plots means the loss of synchronization. 

15 



Consider the Cauchy problem for (j62p with initial condition Ri{x, 0) = ^{x), and the Green 
function G{x, t) such that 

/ + 00 
G{x',t)ip{x-x')dx'. (63) 

-oo 

Let us apply Laplace transform for t and Fourier transform for x 

f'OO f'-\-oo 

G{k,s)= dt dx e-'^+'^''G{x,t). (64) 

Jo J-oo 



G(k, s) = ^ -—. (65) 



Applying (jMj) to (|62|) . we obtain 



Let us first invert the Laplace transform in ()65p . Then, the Fourier transform of the Green 
function: 

/ + 00 
dx e'^'^Gix, t) = e-(2"i+9l'=l")* = e-^'^^'e-^^^^"'. (66) 

-oo 

As the result, we get 

G{x,t) = ((7t)-i/°e-2'^i*L„(x((7t)~^/"). (67) 

where 

1 f + oo 

L^(x) = — dk e-*'=^e-"l'=l" (68) 

27r 7_oo 

is the Levy stable p.d.f. [37]. 

As an example, for a = 1 we have the Cauchy distribution with respect to coordinate 

G(x,t) = - ^f/ ! (69) 

For a = 2, we get the Gauss distribution: 

G{x,t) = (5t)~l/2g-2ait 1 g-xV(4<;t). (70) 

For 1 < a < 2 the function La{x) can be presented as the convergent expansion 

La[x) = > i-x) j sm(n7r/2). (71) 

TTx ^^ n! 

n=l 

The asymptotic (x — )• oo, 1 < a < 2) is given by 

La{x) y (-l)"a;-"" '^^^ V"^ sin(n^/2), x ^ oo, (72) 

TTX ^-^ n! 

n=l 
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with the leading term: La{x) ~ vr ^r(l + a)x " ^, (x — )• oo). 
As the result, the solution of ([58|) is 



Z{x, t) = e*("2*+^(0)) Ly^ + e(<7t)~i/"e-2«i* / La{x'{gt)~^/") if{x - x')dx' + 0{e^) j . 

(73) 
This solution can be considered as a space-time synchronization in the oscillatory medium with 
long-range interaction decreasing as |x|~''"''""'^'. 

For ip{x) = 5{x — xq), solution ([73]) has the form 

Z{x, t) = e^(«2t+e(0)) J^^i/2 ^ £(gt)-i/"e-2"i*L„((2; - XQ){gt)-^''') + 0{e'^)^ , (74) 

and the asymptotic is 

Z{x, t) = e^('^2t+e(o)) j^^i/2 ^ egte-^^^'-K-^Til + a){x - xo)-"~^ + 0(e^)) , x ^ oo. (75) 

This solution shows that the long-wave modes approach the limit cycle exponentially with time. 
For t = l/(2ai), we have the maximum of |Z(j;,t)| with respect to time: 

max |Z(a;, t)| = aj^^ + eg ^^\^ "\ x - xo)'"'^ + 0{e^). (76) 

t>o zTre 

As the result, we have the power law decay with respect to the coordinate for the space structures 
near the limit cycle \Z\ = a^ . 

6 Conclusion 

One-dimensional chain of interacting objects, say oscillators, can be considered as a benchmark 
for numerous applications in physics, chemistry, biology, etc. All considered models were related 
mainly to the oscillating objects with long-range power wise interaction, i.e., with forces propor- 
tional to l/|n — m\^ and 2 < s < 3. A remarkable feature of this interaction is a possibility to 
replace the set of coupled individual oscillator equations into the continuous medium equation 
with fractional space derivative of order a = s—1, where < a < 2, q 7^ 1. Such transformation 
is an approximation and it appears in the infrared limit for wave number /c — )■ 0. This limit 
helps to consider different models and related phenomena in a unified way applying different 
tools of fractional calculus. 
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A nontrivial example of general property of fractional linear equation is its solution with a 
power wise decay along the space coordinate. Prom the physical point of view that means a 
new type of space structures or coherent structures. The scheme of equations with fractional 
derivatives includes either effect of synchronization [8] , breathers |48[ HH] , fractional kinetics [T] , 
and others. 

Discrete breathers are periodic space localized oscillations that arise in discrete and con- 
tinuous nonlinear systems. Their existence was proved in Ref. [51]. Discrete breathers have 
been widely studied in systems with short-range interactions (for a review, see [521 148] ). En- 
ergy and decay properties of discrete breathers in systems with long-range interactions have 
also been studied in the framework of the Klein-Gordon }491 I53j . and the discrete nonlinear 
Schrodinger equations |54] . Therefore, it is interesting to consider breathers solution in systems 
with long-range interactions in infrared approximation. 

We also assume that the suggested replacement of the equations of interacting oscillators by 
the continuous medium equation can be used for improvement of simulations for equations with 
fractional derivatives. 
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